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SUMMARY 


An analysis was made of flight test data to obtain information 
representing the atmospheric turbulence inputs to the test aircraft, The 
basis of the method was the use of the linearized equations of motion of 
the aircraft, with the addition of certain terms representing the aerody- 
namic effects of the turbulence on the aircraft. 

These equations were set up on an analog computer and verified 
by the technique of matching transient response flight data. Then, through 
the use of a feedback system, the desired turbulence quantities were ob- 
tained as system outputs. Inputs to the system were the aircraft response 
and control deflections, obtained from the flight test data. 


The method was found to be valid, assuming the aircraft analogue 


was correct, The resulting data appeared suitable for further analysis To 


determine its statistical qualities. 
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NOMENCLATURE 


Power approach 

Lift coefficient, lift/qs 

Pitching moment coefficient, pitching moment/qS 
Rolling moment coefficient, rolling moment/qSb 
Yawing moment coefficient, yawing moment/qSb 


Lateral force coefficient, lateral force/qS 
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moment of inertia about x axis, slug - ft? 
moment of inertia about y axis, slug - ft? 
moment of inertia about z axis, slug - ft? 
wing area, ft? 

true airspeed, ft/second 

initial airspeed, ft/second 

aircraft weight 


pitching moment of inertia parameter - “у. 
muc 


wing span, ft. 

wing chord, ft. 

mean aerodynamic chord, ft. 
acceleration due to gravity, ft/sec? 
aircraft altitude, ft. 

aircraft mass, slugs 

normal acceleration, az/g 

normal acceleration, ft/sec? 


lateral acceleration, ft/sec? 


= © = rolling velocity, radians/sec. 
= @ = pitching velocity, radians/sec. 
= W = yawing velocity, radians/sec. 


= Time, seconds 
= longitudinal velocity perturbation 
= lateral velocity perturbation 


= normal velocity perturbation 





| INTRODUCT | ON 


In accordance with a proposal submitted to the Bureau of Naval 
Weapons, the James Forrestal Research Center at Princeton University in 
late 1962 undertook an investigation into the lateral flying qualities of 
operational Navy aircraft in the PA configuration approach to carrier 
landings. 

As noted in Reference |, the official proposal, it is typical 
that during this type of approach the lateral stability and control charac- 
Teristics of some high performance aircraft deteriorate with increasing 
angle of attack and decreasing airspeed. This is common in many current 
sweptwing aircraft. The situation is aggravated by the added difficulty 
of approach and landing on aircraft carriers, where the pilot is required 
To maintain a flightpath with extremely narrow tolerances. The difficulty 
is further compounded by natural atmospheric disturbances, and those caused 
by The airflow about the carrier hull and superstructure, which are often 
encountered in the final few seconds prior to touch-down. 

[Т was therefore proposed that a flying qualities research pro- 
gram be conducted at Princeton, utilizing a variable stability aircraft to 
obtain pi lot opinion ratings. The approach to a carrier was To be simula- 
Ted as closely as possible, including effects of the forementioned atmos- 
pheric disturbances. These turbulence effects were to take the form of 
recorded signals applied as inputs to the electronic servos driving the 
control surfaces of the variable stability aircraft. It was therefore con- 


sidered desirable to represent the turbulence as accurately as possible 
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since this factor would probably have an important effect on pilot rating 
in this particular investigation, 

To this end, a suitably instrumented Navy fighter aircraft from 
the Naval Air Test Center, Patuxent River, Maryland, made a series of ac- 
tual carrier approaches under varying conditions of wind velocity across 
The deck and cross wind component. Recordings were made of the aircraft 
response in five degrees of freedom, together with recordings of the var- 
ious control surface deflections. These data were analyzed on an analog 
computer setup of the aircraft equations of motion. By a process to be 
described, it was possible to determine the atmospheric disturbance infor- 
mation in a form suitable for application to the variable stability aircraft, 

This report is concerned with the first phase of this research 
project: the acquisition, reduction, and analysis of the flight test data 


associated with the carrier approaches made by the Naval Air Test Center 


aircraft, 
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[| EQUI PMENT AND PROCEDURE 


|) FLIGHT TEST EQUIPMENT 

The data acquisition flights were flown in an F4B fighter air- 
craft piloted by Navy test pilots from the Naval Air Test Center. The 
F4B (Phantom II) is a two-place, swept-wing, high-performance interceptor 
manufactured by the McDonnell Aircraft Corporation of St. Louis. It is 
powered by two General E!ectric J-79 turbojet engines, each producing a 
maximum of 17,000 pounds thrust at sea level. 

The aircraft was completely instrumented for the data acquisition 
by The Instrumentation Branch of the Flight Test Division, Naval Air Test 
Center. The rear cockpit seat was removed from the aircraft to provide 
space for the instrumentation package. A summary of the flight instrumen- 
tation follows: 

Power Supply: The basic aircraft power package was utilized as 
The instrumentation power source. This consisted of the normal 28 volt 


DC supply, апа 115 volt, 400 cycle AC source. 


Sensing and Transducing: The physical quantities measured by the 


instrumentation system are shown in Table I. Potentiometer type transducers 


were chosen in order to provide compatibility with the telemetering system. 
The potentiometers provided signals of adequate magnitude, and through 
their excitation source also provided zero and full scale reference vol- 
tage for the telemetering unit. 

Table | also lists the basic instrument used in the measurement 
of each quantity. The rate gyros which sensed the various angular velo- 


cities, and the potentiometers which measured the required position angles, 
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were all types commonly used in flight test instrumentation, |t should 
be noted that the actual control surface deflections were measured rather 
than simple stick and pedal deflections, due to the non-linear character- 
istics of the stability augmentation system, Use of the actual deflec- 
tions permitted the aircraft equations of motion to be written in simple 
form. 

Airspeed and altitude were measured by W, O. Leonard force balance 
type transducers. These units were located on the lower shelf of the rear 
cockpit instrumentation sled, as shown in Reference 2. The transducers, 
operating on the force balance principle, provided linear outputs for air- 
speed and altitude without the use of cams, function gears, shunts, efc. 
The force sensors were isolated from the output devices, thus minimizing 
response to transients. As a result, the units were characterized by a 
pressure accuracy ot up to .1%. In the case of the altimeter this meant 
accuracy to within three inches for a two hundred foot aircraft altitude. 

Lateral and normal acceleration were sensed by Donner model 4310 
linear accelerometers installed with a range of + .5 G and + 2.0 G, respec- 
tively. These units exhibited natural frequencies much higher Than The 
Test aircraft, thus precluding resonance problems, Their response was es- 
sentially flat to 8 cps. It was determined, however, during installation, 
that some H. Р. system noise was being passed by The acceleromeTer。 In 
order to eliminate this noise, a low pass filter was therefore designed 
and installed between the accelerometer output and the telemetering unit. 

The filter, designed for a 5 cps corner frequency, eliminated 


The high frequency noise but provided essentially flat response over the 
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expected aircraft response range, which was estimated to be no higher than 
l cps, well above the Dutch roll frequency. 
A schematic of the filter circuit is shown below, together with 


a sketch of the filter frequency response characteristic: 
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Telemetering and Recording: Telemetering equipment of the pulse 


width type was installed in the aircraft as part of the instrumentation 
package. This equipment, manufactured by the Applied Science Corporation 
of Princeton (ASCOP), was the manufacturer's D Series Pulse Width Multi- 


coder-Telemeter, capable of sampling up to 43 channels at a sampling rate 





of 20 times per channel per second, using a Model DSG20 Commutator-Gear 
Box. 

Since each channel was sampled 20 times per second, it was felt 
that good fidelity would be obtained with input frequencies as high as 
the 2 to 5 cps range. Many of the measured quantities were sampled on each 
of two or three channels, and in these cases it was felt that good repro- 
duction was assured at all frequencies encountered in the investigation. 
Table 11 lists the measured quantities and the number of sampling channels 
used for each. For those quantities sampled in more than one channel, the 
sampling was, as far as possible, equally spaced in time. 

The sampled data were then converted into PW data by The Model 
DKH-! Keyer. 

The pulse width information was then recorded directly in the 
aircraft by means of an Ampex Model 800 tape transport which was mounted, 
together with its associated power supply, in the rear cockpit installation. 

Provisions were also made for simultaneous voice recording by 


the test pilot. 


2) GROUND EQUIPMENT 

The following units were available at the Forrestal Research 
Center for analysis of the flight test data: 

Magnetic Tape Recorder: An Ampex Model 309C dual track, three 
speed tape recorder was available for playback of the recorded flight data, 
and also for use in the preliminary analysis studies which were conducted. 


Telemetering Ground Station: The telemetering unit used for de- 
coding the taped pulse width flight test data was the ASCOP M Series PW 





Ground Station. This unit continuously sampled the recorded pulse width 
information and converted this data into continuous voltage form. The 
information in this form was then compatible with the analog computer for 
analysis purposes. 

GEDA Analog Computer: The basic tool in the analysis process 
was the Model L5 (GEDA) linear electronic differential analyzer, manufac- 
tured by the Goodyear Aircraft Corporation. This computer provi ded twenty- 
four automatically stabilized DC computing amplifiers, each with an open 
loop gain of greater than 5 x 10”, and with negligible drift. The computer 
incorporated an automatic error indicator, and had a guaranteed accuracy 
of I$. There were provisions for accurately setting computer board po- 
tentiometers using a special null indicator. 


Data Recorder: Visual presentation of data was made possible 





by a Sanborn Model 154-1008 four channel recorder. This instrument had 
a very high natural frequency (42 cps), and low frequency response flat 


to zero cps. 


3) INSTRUMENT CALIBRATION 


Calibration of the flight test instrumentation was conducted by 
the {Instrumentation Branch of the Flight Test Division, Naval Air Test 
Center. The calibration charts are shown in Reference 2. All gradients 
were determined to be essentially linear within the range of measured data 


encountered during the test flights. 
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4) PROCEDURE 


Data Acquisition: 
in December, 1962, the instrumented F4B aircraft, piloted by a 


Navy test pilot, made a series of thirty-three landing approaches to the 
aircraft carrier USS Enterprise. These approaches and landings were made 
over a span of five separate flights. The heading and speed of the car- 
rier were varied during the course of the series to provide a number of 
different conditions of wind velocity across the deck, and cross wind com- 
ponent. The wind conditions for each approach are shown in Table lil. 

ІТ was expected that the most severe turbulence conditions would 
be experienced in the presence of a "right to left" crosswind. This was 
due to the extremely turbulent flow about the carrier "island" being direc- 
ted into the path of an approaching aircraft, usually at the critical point 
occurring a few seconds prior to touchdown. It was felt that by making 
the approaches under varying wind conditions, it would be possible to ob- 
Tain turbulence data which would be representative of the full range of 
severity normally encountered. As noted in Reference 3, the turbulence 
"burble" associated with the Enterprise was probably the extreme, due to 
The very boxlike design of the ship!s "island", 

The Naval Air Test Center provided facilities aboard the Enter- 
prise for tape playback after each test flight. On the spot assessment of 
data quality was thus possible, and minor adjustments were made, as neces- 


загу, in the instrumentation, 
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Aircraft Calibration: 


Since the basis of a successful analysis of the flight test data 
was to be the setup of the aircraft equations of motion on the analog com- 
puter, it was necessary, as the next step, to obtain a valid set of equa- 
tions representing the F4B aircraft in the PA configuration. The basic 
method used was the matching, on the computer, of transient response flight 
test data. The method will be discussed in some detail in a subsequent 
section. 

To this end, two additional flights were made of the test air- 
craft for the purpose of obtaining transient response data. The aircraft 
was flown in the PA configuration, and pilot control inputs of the pulse, 
step, and doublet types were made, The aircraft response was recorded 
as before. 

In order to prevent stability augmentation equipment inputs to 
the system, these units were made inoperative during the transient res- 
ponse data runs, 

Data Analysis: 

Having thus acquired the necessary flight test data, the actual 

analysis was then conducted at the Forrestal Research Center, This analysis 


is discussed in detail in the following section of this report. 
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[|| | DATA ANALYSIS: DESCRIPTION AND DISCUSSION 


The first step in the data analysis was the determination of a 
valid set of differential equations representing the F4B aircraft in the 
PA configuration. It was assumed that the motions of the aircraft while 
approaching a carrier in this configuration could be described by the usual 
linear equations. This assumption seemed valid since the final part of 
a typical carrier approach, where the turbulence effects are normally most 
severe, is basically a constant power, constant airspeed, shallow descent 
to a touchdown; the aircraft motions are approximately perturbations from 
the steady state glide path, 

The resulting equations are shown below, written in dimensional 
form for compatibility with the flight data and analog computer. For con- 
venience, the terms of each equation have been divided by the coefficient 
of the highest appropriate derivative. The assumption was also made that 
the aircraft forward velocity remained essentially constant during the 
runs, thus eliminating the drag equation, and also eliminating longitudinal 


velocity terms from the remaining equations. 


Longi tudi nal: 


m C 
(1) lift eguation: a-0+[ Cra ]a +[ “551 65 = 0 
20-7 2T 


(2) pitch equation: 0 - [ Став ) = [ Ed 16 = [ Cm. 16 
HT HT HT? 


- [ ms, 10 £ D 
H T? 


















a ee A a oc NN у‏ كت 


n ID L! Ii 


. í в 


o ж - ж жоң кк ж. " [ie m ' 


に -- ወ ይ. መሽ お = > | am + ) -! - - 


m E A жн кн д 
-— 020 - ም اسه > ا‎ шш — ብ 
A. SI шш - a -— ቐ 8117 T 


ай а. >» 1 س‎ à PAX o ән 


|o — We 

N — £n 
TB. FAA . US C ፦. 
Gail; Ó, w < | p 
ーー ニー ザー не ӨТІ «аі аса A 
БАЛАА ——— |) 
CO  . 14. (^ m че э eem 
Ааа se 


fm! aa Gan 









EIS 


(5) normal acceleration: а = -C, ) Vo Ja + С, ( 6ك‎ 
DT 5 2T 


Lateral: 
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These uncoupled sets of equations were then separately mechan- 
ized on the GEDA analog computer. Schematic circuit diagrams of the two 
setups are shown in Figures | and 2. 

All factors in these equations were known accurately, except 
for the various stability derivatives. Approximate values for the deriva- 
tives, for the PA configuration, were obtained from a report of the 
McDonnell Aircraft Corporation. These values determined the initial po- 
tentiometer settings for the computer setups. Shown below is a list of 
the various quantities utilized in the initial formulation of tne aircraft 


equations for the PA configuration: 
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The equations of motion based on these values are as follows: 
Longitudinal: 
- Ó + ,406 a + .0456 8, = 0 


Qe 


8 + ,162 8 + ,067 ሴ + .357 a + .792 6, = 0 


n+2.81 a - .311 55 =0 


Lateral: 
- B - .0832 В - .990 ў з .0109 ф -.Із9 ф з .0166 6. = 0 
-9-1.25$ *.986 q - 12.2 B - .510 6, + 2.80 ба > 0 
= 9 


- 4 - .264 $ + .0288 ġ + 1.89 B - .693 ô. - .0556 ôa = 


Improving the validity of these equations as representations of 


the aircraft thus requirod a more accurate determination of the stability 
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derivatives. This was accomplished by a transient response curve fitting 
technique, such as outlined in Reference 4, and described in detail in 
Reference 5. The mathematical basis for the validity of the method is 
found in Reference 6, 

For this method a number of the transient response runs were 
selected by reference to Sanborn recordings of the aircraft response. The 
selection was made on the basis of freedom from telemetering dropouts, and 
relative freedom from apparent extraneous inputs, such as turbulence effects, 
These response curves were of the rudder pulse, rudder doublet, aileron 
step, aileron doublet, stabilator pulse, and stabilator doublet types. 

Transparent overlays were made from the Sanborn records, incor- 
porating the scale changes necessary to provide agreement with the scaling 
of the computer outputs. It was then possible, by means of the ground 
station equipment, to apply voltages to the aircraft analogue represent- 
ing the actual aircraft control motions, as forcing functions to either the 
longitudinal or lateral sets, as appropriate. Sanborn recordings of the 
analog responses were made, all to the same scale as the transparent over- 
lays. The two sets of data were compared, and by an iteration procedure, 
the potentiometer settings representing the various stability derivatives 
were changed until close agreement was attained between flight response 
and analog response. 

In the lateral case the matched quantities were rol! rate and 
yaw rate. It was found that a close match to the flight data could be 
fairly rapidly obtained by corrections to those derivatives which had major 


etfects on the important lateral response parameters, i.e., Dutch roll 





Sube 


period and damping, initial roll and yaw magnitudes, spiral stability. 
Improvements To the initial match were obtained by a systematic iteration 
procedure, varying the potentiometer settings in turn until the desired 
match resulted. 

In the initial study of the lateral responses it was immediately 
apparent that the flight data was not ideal. During even the best runs, 
the aircraft obviously had not been in steady flight at the instant of the 
control inputs; there were initial conditions which had to be consi dered 
in the matching process. By means of the analog initial condition cir- 
cuitry it was possible to apply to the computer, as initial conditions, 
the actual flight data response in roll rate and yaw rate existing at the 
instant of control deflection. This was considered to be an improvement 
but it was noted that since the lateral equations represented a fourth 
order system, two additional initial conditions, B and q, for example, 
were necessary for accurate solution, These were not available and some 
error was thus to be expected. In the iteration procedure it was found 
that the best matches were obtained with flight data where the aircraft 
was very close to steady state flight conditions at the time of control 
deflection. This was to be reasonably expected since in these cases all 
initial conditions were approximately zero. 

1ተ was also apparent from observation of the flight data that 
some atmospheric turbulence existed during the tests, This, of course, 
could not be accounted for in the matching process, and was thus a source 
of error. For this reason the smoothest runs were selected for the match- 


ing process, 
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It was noted also that during the various transient response runs 
there were inadvertent inputs of those control surfaces which theoretically 
were being held in the neutral position. These inputs were also played into 
the analog. 

ІТ сап Бе seen from Table | that only the left hand aileron and 
spoiler were instrumented for deflection measurement. Since the equations 
of motion as formulated on the computer were written in terms of aileron 
inputs, which are the sum of right and left aileron deflection, it was 
necessary to derive some function of left spoiler deflection which might 
be utilized as equivalent right aileron deflection. To this end, a cali- 
bration run was made to determine the relationship between aileron and spoiler 
deflection. The results were of the form shown in Figure 3. The relation- 
ship between the two was shown to be approximately linear. Therefore, it 
was possible to express right aileron as a linear function of the recorded 
spoiler deflection. The aileron signal fed to the analog was: 

(7) 6. = ба, + Š = ба, Za Š sp 


On L 


IT can be seen from Figure 3 that each aileron exhibited a small 
upward deflection range of approximately three degrees. This property was 
ignored in determining the proportionality constant between aileron and 
spoiler. |t was felt that the upward deflections were primarily for the 
purpose of assisting the effect of the spoiler deflections on the same 
wing. Іп the slightly open position the spoilers were assumed to be rela- 
tively ineffective in disrupting the airflow, and the small upward aileron 


deflection made up for the deficiency, thus maintaining system linearity. 
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By means of the analog limiter circuits, the voltages representing these 
upward aileron deflections were eliminated from the computer inputs. 

[T should be noted that in spite of the turbulence and initial 
condition problems, it was possible to obtain good matches between flight 
data and analog output. The accurate determination of all the various 
stability derivatives would probably have been difficult because of these 
problems, but the primary interest here was the matching of aircraft re- 
sponse, This could be done by matching the important lateral parameters 
mentioned previously. Since these parameters are primarily controlled by 
relatively few of the derivatives, it was possible to obtain a good match 
even Though some of the derivatives might have been in error. 

It should also be noted that since the instrumentation used in 
obtaining data for this matching was the same as subsequently used during 
the actual carrier approaches, then to the extent that any instrumentation 
errors would have caused false matching, the effect would tend to be can- 
celled during the actual approach data acquisition. 

Shown below are the circuits utilized to provide inputs to the 
analog from the recorded data. By means of appropriate potentiometer set- 
tings the quantities were all scaled correctly for compatibility with the 
analog setup. Biasing of all quantities to the correct zero position was 
done by adjusting the biasing knobs of the ground station while referring 
fo voltmeter readings of the various tape outputs. 

Shown in Figures 4, 5, and 6 are some of the matched transient 
response curves. It is to be noted that the aileron step inputs were used 


primarily to match initial roll rate amplitude. It was felt that the aircraft 
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was out of the linear range shortly after application of an aileron step 


and this precluded any additional useful matching. 
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The longitudinal responses were matched by comparing pitch rate 
and normal acceleration. The method utilized and the problems encountered 
were similar to those discussed in reference to the lateral equations. 

The phugoid mode was not visibly excited in any of the transient runs and, 
Therefore, the matching was confined solely to the short period oscil- 
lations, The matching iteration procedure converged much more rapidly than 
in the lateral case. This was to be expected since in this longitudinal 
approximation there were only two degrees of freedom and, consequently, 


fewer controlling derivatives to be adjusted. 
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shown below are the input circuits to the longitudinal analog: 
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In Figures 7, 8, and 9 are shown some of the matched transient 
response curves. 

It should be noted that the initial matching, both longitudinal ly 
and laterally, showed a number of discrepancies between the actual aircraft 
response and the analog response which was based on the stability deriva- 
tives obtained from the McDonnell! report. Іп the lateral case the analog 
response showed a predicted Dutch rol! period of approximately 4.5 seconds, 
with neutral damping. As can be seen in Figures 4 and 5, the aircraft 
actually exhibited a well damped Dutch roll of about a 3 second period. 
Longitudinally, the major discrepancy was in the apparent degree of air- 
craft static stability, which was considerably greater than the predicted 
value. 

These discrepancies are illustrated more clearly by reference to 
the final equations of motion which resulted from the matching process. !f 


these equations, which are shown below, are compared term by term with the 
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original set, the changes which had to be made during the matching are 


apparent. 
Longitudinal: 
(8) 4-0 +.472 4 + .0436 6, = 0 


+ „858 6 + .067 а + 9.00 с + 4.50 6g = 0 


Ds 


(9) 
MN + 3.25 a - .311 6, = 0 
Lateral: 
(11) - B - .0956 B - .990 $ + .0110ф + .140 ф + .0209 6, = О 
(12) - ዌ - 2.48 ቆዕ - .984 y - 15.8 B + .497 5, + 5,48 57 =0 
№5) - y = ጋጋሪ ў + .0288 о + 3.38 B - .674 6r - .061 Sa = 0 
|t can be seen that in the lateral set the directional stability 
and yaw damping, which are the primary parameters affecting Dutch roll 
period and damping, respectively, were each approximately doubled in the 
matching process. M Longitudinally, it can be seen that the static stability 
was also increased considerably during the matching. Although the high 
value of Сто, which finally produced a good match, Was probably unrealistic 
it at least indicated that the predicted value was too low. It is noted 
again that accurate determination of each stability derivative was by no 
means assured. However, the aircraft response curves were fairly accurately 
matched, and this was the important point for purposes of this investigation. 
With the equations of motion of the aircraft thus formulated, 
the next step was the actual analysis of the flight data to obtain the 
turbulence information, Preliminary studies were made to determine the 
feasibility of the proposed method of solution. For these studies the 


lateral equations of the Navion aircraft, as obtained from Reference 5, 





-2l- 


were set up on the analog computer. The basis of the method was the well- 


known feedback principle: 





From the above diagram, C/R з А/(І - Ав). If A, the forward 
loop gain, is very much greater than unity, then, 
C = 1 
R B 
If B were an aircraft transfer function, e.g., q/6a, then, 


Cw LL =< да 


R Ф/ба Ф 
ІТ seemed possible, then, that if known aircraft response in roll rate, for 
example, were fed into such a circuit at "R", where "B" was the aircraft 
fransfer function as set up on the analog computer, then the circuit out- 
put at "C" would be a voltage representing the ai leron deflection which 
had caused the roll response. 

The validity of this theory was tested by recording on tape the 
analog response in Ф, ў, and B, to an arbitrary sequence of voltages repre- 
senting aileron inputs to the system. This was accompli shed by connect- 
ing a potentiometer to the appropriate points of the analog circuit and, 
in effect, "flying" the aircraft in the lateral modes. The previously 


discussed feedback circuit was then connected as shown below: 
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The taped response in roll rate was then played into the circuit at "К" 

and the output "C" was compared with the original aileron signals. The 
match was excellent with forward loop gains as low as 10, Below this value 
the reproduction decreased in accuracy since (C/R) = (1/B) was no longer 

a good approximation, 

The next logical step was to study the effects of two control 
inputs, both aileron and rudder. Again the analog was "flown", This Time 
with aileron and rudder. The same aircraft response quantities were again 
taped. The analysis circuit was similar to the previous case, with the 
addition of another loop. The circuit is shown below in block diagram 
form. The taped responses in Фф апа y were played into the circuit at R 
and R' respectively, and the outputs at C and C' were compared with The 
original 6, and 6-. Again the match was excellent, down to forward loop 


gains of IO. 
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Figure IO shows a more detailed schematic of this last analysis 
circuit. |t was important to provide negative feedback in the loops in 
order to maintain system stability. 

The results of this preliminary study pointed up the uniqueness 
of the solutions to the aircraft equations of motion. A given combination 
of aileron and rudder inputs results in a unique combination of roll, yaw, 
and sideslip response. 

Originally, it had been proposed that this exact analysis method 
be employed in the study of the flight test data. In the lateral case, 
for example, aircraft response in roll rate and yaw rate were to be applied 
to a similar analysis circuit, where the analog would be of the F4B. From 
the resulting outputs of 6, and Š, Then would have been subtracted the 
recorded flight test control deflections. The difference would be an 


aileron and rudder signal representing the turbulence input to ተከፀ የ48, 
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This information, suitably modified to account for the difference in con- 
trol power between F4B and Navion, might then be applied to the variable 
stability aircraft as a simulation of the turbulence input. 

The longi tudinal analysis would be similar. Aircraft response 
in pitch rate would be applied to the appropriate circuit and the output 
65 would yield an equivalent ôs representing the turbulence input. 

This method, however, obviously neglects the linear perturbations 
of aircraft motion due to turbulence, those motions along the aircraft 
axes. lf assumes that the atmospheric turbulence effects may be des- 
cribed accurately by the resulting moments about the aircraft axes. For 
purposes of this research project reasonably valid simulation may be at- 
tained based on this concept. 

However, in the interest of more accurate gust simulation, it was 
decided to conduct some initial research into the possibility of obtaining 
turbulence data which might be expressed in five degrees of freedom, thus 
including any lateral and normal acceleration effects, In addition, the 
possibility existed that some qualitative information concerning the char- 
acter of the turbulence itself might result from such a study. 

It was known that it is possible to visualize the atmosphere 
Through which the aircraft moves as being in motion at some steady mean 
value, with random turbulent fluctuations superimposed. These random velo- 
cities may be considered as inputs to the aircraft linear system, which 
may respond in six degrees of freedom. The velocity inputs are assumed 


small, permitting linearization of the resulting aerodynamic forces. 
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As outlined in Reference 7, it is convenient to introduce these 
inputs into the equations of motion in the form of velocity components 


ug, Vo and Wg (subscript "g" standing for gust). If it is assumed That 


g 
the gust wavelengths are all large relative to the aircraft dimensions, 
Then The gradient, or velocity variation, of these gusts across the length 
and span of the aircraft may be considered essentially linear. IT can 

be shown that under this assumption the effects of the gust components and 
their gradients are aerodynamically equivalent to perturbations in angle 
of attack, sideslip, roll rate, yaw rate, and pitch rate. These are noted 
аз ас, Ва? Pgs Те Og. 

The aircraft equations of motion incorporating the modified var- 


tables thus become: 


Longi tudi nal: 


° : C 
(14)9 70 + e EDE ts ôs = 0 
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With the gust terms rearranged to point out their role as forcing functions, 


the equations are written as: 


Longi tudinal: 
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It is to be noted that these equations are not valid for short 
wavelength gust components. 

However, it can be shown that the effect on aircraft response of 
the higher frequency gust components is negligible. For example, the ef- 
fective gust frequency, Og, may be expressed by: 


UR‏ دن 
di w‏ 


where à is the gust wavelength, For a short à, say equivalent to the air- 
craft wingspan, Wg is on the order of 5 cps, for the F4B in the PA confi gur 
ration. This is very much higher than the aircraft natural oscillatory 
frequencies and consequently the aircraft response is negligible. 

In addition, for gust wavelengths of this order, and shorter, 
the aircraft begins to "average out" the gust velocities, with the result- 
ing net gust acting on the airframe being near zero. Thus, for these higher 
frequency components there is a compounded attenuation - that due to the 


aircraft dynamic characteristics, and that due to The cancelling effect 


AL] 


IP a 


just described. 
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ІТ should be noted that if the F4B were to respond dynamically 
to the higher frequencies, this last effect, that of cancelling, would 
cause some error if the acquired turbulence information were applied to the 
study of an aircraft differing greatly in size from the F4B. It is quite 
possible that the aircraft in question might respond to those gusts which 
were cancelled by the F4B and were thus missing from the acquired data. 
However, all aircraft which will be simulated by the variable stability 
Navion, for this investigation, are similar in size to the F4B, and there- 
fore could be expected to "average" the gusts in the same manner. 

With these modified equations of motion in mind, it was desired 
to obtain, as a solution, the various gust quantities. The method, although 
based on the previously discussed feedback principle, can probably be best 
described by a single degree of freedom example. For this purpose it is 
assumed that the aircraft can respond only in roll, and that the only gust 


quantity present is фа. The equation of this system is then: 





(24) Фф - “р о ë Wels, Ds us 64s ôr -~ こ %p op。 = 0 

23,7 Jx T° ጀው 2109 E 

This equation is set up on the analog computer and a feedback 
circuit is connected as shown on page 29. 

The output of the forward loop amplifier is defined as - фа, 
and this quantity, modified by the appropriate proportionality constant, 
is used as an analog forcing function. The recorded values of ~ ôr and 
- 5, are also used to force the computer setup. The following equation, 


which expresses the aircraft response, is satisfied at point A; 


recorded = Pcontrol response * Pgust response 
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Assuming that the circuit is dynamically stable, then the - фа 
quantity actually appears at the amplifier output. This circuitry is simi - 
lar to those previously di scussed. 

The more complicated circuit which was proposed to solve the en- 
Tire lateral set is shown in detail in Figure Il. It is, however, based 
on the same ideas. The three equations which are satisfied at the forward 
loop amplifier inputs are: 


B recorded ч Ва 7 Boust response ` Pcontrol response ^ 


Vrecorded ^ $gust response ` $control response ^ 0 


Vrecorded ” Waus+ response ` "control response ” 0 

It is to bé noted that in the case of B, the recorded quantity 
from the tape was based on an accelerometer measurement and, thus, con- 
tained the quantity Bg, in addition to Baus+ response? the aircraft re- 
sponse to Ва. Therefore, to satisfy the above equation, the Bg was fed 


back to the amplifier input. In the case of d$ and |, the aircraft response 


only was included in the measured response since the sensors were rate 





250% 


gyros. Approximately infinite resistance feedback was therefore provi ded 
around the forward loop amplifiers in these two channels. 

lt was di scovered, however, that the circuit shown in Figure 
Il was not satisfactory, During operation of the circuit the output quan- 
tity - Bg exhibited a divergence. This indicated some positive feedback 
within the multi loop system, It can be seen that there were nine loops 
through the aircraft analog. Г seemed probable that the - و8‎ output was 
controlled primarily by the loop through amplifier 9 which is shown below 


schematically: 


| rcarded 





This circuit obviously provides positive feedback and is Thus 
unstable. This problem emphasized the point that although the equations 


were physically satistied in the circuit as shown, it was necessary To 


insure system stability in the electronic sense. A study of the circuit 


showed that there were additional unstable loops, but that their effect 
was minor compared to the effect of the stable loops which fed back To 
the same points. 

In order to provide stability, the analysis circuit was modi fied 


The difference here is that the - Bg quan- 
The 


To the form shown in Figure 12. 
tity input to the analog is now taken directly from the recorded data. 


analog circuit was modified to maintain satisfaction of The equations of 
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motion. Circuit stability was now attained since the unstable - Ва loop 
through amplifier 9 was now eliminated. 

To test the validity of the analyzer circuit, a preliminary study 
was made in which the aircraft analogue was forced by a ficticious Bg quan- 
tity as obtained from a potentiometer output. The aircraft response in 
roll rate, yaw rate, and sideslip was recorded on tape. It is to be noted 
that the recorded sideslip was made the sum of the analog B response and 
the forcing function, Bg. This was to provide similarity to the actual 
flight data where, as previously mentioned, the B response was based on 
an accelerometer measurement, These recorded responses were then played 
into the analyzer circuit in an attempt to duplicate the original forcing 
function, 

Figure 13 shows the original forcing function and the recorded 
aircraft responses. When these responses were played back into the апа- 
lyzer circuit, the Pg and Wo outputs were essentially zero, as expected. 
The Ва output was found to be a good reproduction of the input wave form, 
although there was a very low frequency oscillation superimposed on The 
trace, This is shown in Figure 14(a), where the solid line is a tracing 
of the original Ва input. It was determined experimentally that this low 
frequency drift was strongly affected by the biasing adjustment of The 
tape output from the ground station equipment, Careful adjustment of the 
bias practically eliminated the drift; but since this method would not be 
feasible with actual flight test data, another method was sought. IT was 
decided to utilize a high pass filter to eliminate the low frequency о5- 


cillation. This was connected as shown: 
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Since the use of such a filter would result also in the attenua- 
tion of some of the higher frequencies, the corner frequency was selected 
on the basis of compromise between minimum attenuation and minimum diver- 
gence. The results are shown in Figure 14 (b), (c), and (d). A corner 
frequency of 1/44 cps, as provided by R4 = 7 በ69 ©, seemed to give the 
best reproduction, 

Throughout this discussion it has been assumed that the aircraft 
response data were available in the form of В, Ф, and ў. IT can be seen 
from Table | that lateral acceleration, ay, was recorded during the flight 
Tests, rather than sideslip. It was necessary to obtain B from this data 


by means of the following form of the side force equation: 





(25) + B = m لك‎ aj = Сур b の - Cy. b i - Yo ô 
ረመ ЕЕ о 7 


The summing and scaling circuit utilized to solve this equation 


for application to the analyzer is shown below: 
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Since an accurate representation of B Та алан on reliable values 
ОЁ Сув, Су. Cys and Cys _, and since these derivatives were not involved 
in the previous matching process, it was necessary to check their validity. 
The method used was to compare the quantity ay from the flight data with 
the derived ay obtained by summing the appropriate analog outputs in re- 
sponse to a rudder pulse. The results are shown in Figures 4, 5, and 6. 
It was felt that the matching was accurate enough to justify the use of 
B, as obtained from equation 25. The resulting equation for ay is: 


(26) + بت‎ - 25 Сув В +.00/2)У25 b Cyp ፍቶ ፈ፻ 2)Vo $ b Cy ў 
т т 2Vo m ZVo 


+ (0/2)%o SC). бе 
m 
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The summing and scaling circuit for obtaining ау for the matching 


process is shown below: 


2 





» HA 


The analyzer circuit designed for use with the longitudinal equa- 
tions was based on the same principles which governed the lateral circuit. 
The schematic of this circuit is shown in Figure 15, The following equa- 


tions were satisfied at the inputs to the forward loop amplifiers: 


° e 


e 
= gust response * ®control response 


aircraft response 


| - + 
“aircraft response “ “gust + %gust response * “control response 


It should be noted that the reliability of the analyzer output 
was probably improved if the aircraft maintained a very steady, smooth 
flight path during the carrier approaches. The reason for this is that 
with the aircraft responses very small, the aircraft analogue, in The ana- 


lyzer feedback toops, was relatively inactive. This was beneficial since 





the analogue had known inaccuracies. |t was therefore desirable that the 
aircraft stability augmentation equipment be operative as an aid in main- 


taining steady flight during the carrier approaches. 
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IV RESULTS 


The next step in the data analysis was to play into the analyzer 
circuits the appropriate quantities from the flight test tape recordings. 
Of the five flight test tapes, one was determined to be virtually unusable 
due to superimposed noise. The four remaining tapes were each played into 
the analyzer and Sanborn recordings were made of the five gust quantities 
for each carrier approach. Also recorded for reference purposes were the 
aircraft response in five degrees of freedom and the various control 
deflections. 

Figures 16 through 31 are reproductions of gust quantities from 
the lateral and longitudinal analyzers. Also shown are traces of aircraff 
response and control deflection. Aircraft altitude also appears on the re- 
cordings as an aid in orientation. It was apparent from the traces that 
some drift was still present in the Bg channel, probably due to incorrect 
biasing at the ground station output, Further increase of the filter cor- 
ner frequency was not practical due to the extensive attenuation of The 
higher frequency components. 

Those approaches which are missing from the set of Sanborn re- 
cordings Were unusable due to the effects of "dropouts" during the analyzer 
operation. As mentioned previously, the tape of the second flight, con- 
taining approaches 8 through 13, was completely unusable due to poor quality 
of recorded information, 

Study of the traces revealed a number of points, IT was expec- 


ted that if the analyzer circuits were completely valid the output gust 
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quantities would be completely random in character, reflecting the nature 
of the actual turbulence. However, Figure 19, for example, shows that 
there were occasional sections of the various gust traces which were perio- 
dic in character. Conversely, many of the remaining traces, Figure l6 for 
example, appear quite random. Some of the oscillations appearing were of 
the approximate character of the aircraft short period and Dutch roll oc- 
cillation. The implication here was that the F4B analogue was not com- 
pletely accurate. As a result, some of the aircraft response characteris- 
tics appeared in the analyzer output. 

In spite of the apparent imperfections in the turbulence data, 
some qualitative information was immediately available. The records showed 
that for all approaches there was a marked increase in the magnitude of The 
gust quantities in the few seconds prior to aircraft touchdown. This was 
a verification of the existence of the turbulent "burble" immediately astern 
of an aircraft carrier deck. 

In order to detect any additional trends in the turbulence quan- 
itites, the RMS value of each quantity for each approach was determined 


with the aid of the following circuit, which included the Donner Model 


3751 Function Generator: 
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t 
(27) RMS = | X p Kis 


In order to eliminate the effect of any DC bias or low frequency 
oscillations, the various gust quantities were channeled through a high 
pass filter prior to entry into the above circuit. A corner frequency of 
.20 cps was selected for the filter, since it was felt that frequencies 
below this would probably not have a significant effect on pi lot opinion. 

The RMS values of the various quantities were tabulated and are 
shown in Table Il. Omissions in the Table were due to the effects of tele- 
metering dropouts, These tabulated values do not include the effect of the 
final "burble", No apparent trends were determined from study of the RMS 
tabulation, except for an indication of the general level of each quantity. 
These levels remained approximately constant during the course of the ap- 
proaches, indicating that atmospheric conditions also generally remained 
the same. 

It was also noted, with reference to Table 111 that the angular 
velocity gust quantities, Фа» fo; and з were all of approximately the 
same magnitude. The angular quantities Bg and aq were also approximately 
equal. This was a strong indication that the turbulence away from the 
imnediate vicinity of the ship was indeed isotropic, as is normally assumed. 

It was expected that the magnitude of the gust quantities within 
the "burble" would be a reflection of the existing relative wind conditions 
over the carrier deck, Although no formal tabulation was made of the RMS 


values of the gust quantities within the "burble", study of the traces 
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revealed that these RMS values were approximately doubled in the final 
three seconds of each approach. There was no apparent difference in the 
character of the "burbles" with varying relative wind condi tions. 

[ተ was noted, however, from study of the gust traces that the 
"burble" was initially encountered at varying distances from the ship. 
Figure 32 shows the point of entering the "burble" ftor all of the approaches 
represented by Figures 16 through 31. |t is interesting to note that those 
points furthest from the carrier were in the presence of high relative 
wind, with no x-wind component, Those points nearer the ship were with 
less wind velocity. Finally, those points appearing 2 and 3 seconds prior 
to touchdown were in the presence of "right to left" x-wind, indicating 
that the turbulence trail was being blown across the flightpath at This 
point, 

In Figure 33 are plotted the pertinent lateral mode data from a 
typical approach. The increase in the magnitude of the gust quantities in 
the final few seconds before touchdown is apparent. However, it was also 
noted that the gust quantities at this point had some periodic characteris- 
tics. |n addition, the oscillations are approximately of the same period 
as the aircraft roll and yaw response. In fact, in the case of ф апа Фа» 
the oscillations are exactly 180” out of phase. Because of apparent "im- 
purities" such as this in the turbulence data it was felt that a more de- 
tailed statistical analysis of the data should be made to determine ¡Ts 
true spectral character. Such an analysis is beyond the present scope 


of this project. 
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V CONCLUSIONS AND RECOMMENDAT I ONS 


Based on the results of this investigation, it is concluded that: 
|. Accurate representations of the atmospheric turbulence may be obtained 
by the analysis method described in the report, assuming that the aircraft 
analogue is valid, 

2. The consistent similarity in the magnitude of the angular velocity 
gust quantities, and the similarity in the magnitudes of gust angular quan- 
tities, as shown in Table |V, indicated that the turbuience away from the 
immediate vicinity of the carrier was isotropic in character. 

5. The atmospheric turbulence in the vicinity of an aircraft carrier is 
characterized by a strong "burble" in the vicinity of the approach end of 
the carrier deck. It is not possible, based on these results, to give 

a quantitative description of this turbulence although it appeared that The 
character of the "burble"did not vary noticeably with changing relative 
wind conditions, However, the point at which the "burble" is initially 
encountered is a function of the relative wind conditions, 

It is recommended that: 

l. A detailed spectrum analysis be made of the turbulence quantities re- 
sulting from this investigation, in order that undesirable portions of Тһе 
spectrum be eliminated, making the information suitable for application 

to the variable stability Navion. 

2. A further statistical analysis of this data be made to determine quan- 
titatively the character of its components, 

5. An investigation be conducted into the applicability of this analysis cir- 


cuit for the general solution of complex inverse transfer function problems, 
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TABLE | 


RECORDED FLIGHT TEST DATA QUANTITIES 


QUANTITY 

Roll rate, c.g. 
Pitch rate, c.g. 
Yaw rate, c.g. 


Roll rate, seat 


Pitch rate, seat 


Yaw rate, seat 


Pitch and roll rate, 


vertical gyro 


Lateral acceleration 


Normal acceleration 


Rudder deflection 


Left aileron defl. 


Left spoiler defl. 


Stabi lator defi. 


Rudder pedal defi. 


Lateral stick defl. 





az 


бг 


Šs 


бред 


Ostick, lat 


SENSOR. 
Humphrey rate gyro, Model 01-0258 
Humphrey rate gyro, Model 01-0238 
Humphrey rate gyro, Model 01-0238 


Minn.-Honeywell rate gyro, Model 
JG 7005A4-9 


Minn.-Honeywell rate gyro, Model 
JG 7005A4-9 


Minn.-Honeywell rate gyro, Model 
JG 7005A4-9 


Mi nn. -Honeywell rate gyro, Model 
JG 7044A4 


Donner linear accelerometer, 
Model 4310 


Donner linear accelerometer, 
Model 4310 


2 KQ, modified, spring loaded, 
selsyn potentiometer 


IRC 2 KQ dual, spring loaded, 
potentiometer 


Bournes 4 inch linear 
potentiometer 


2 K Q, modified, spring loaded, 
selsyn potentiometer 


spring loaded, modified selsyn 
potentiometer 


spring loaded, modi fied selsyn 
potentiometer 


. wor ТТ жен», 


мы 





TABLE | 
(Cont i nued) 


UANTITY SYMBOL SENSOR 


spring loaded, modified selsyn 


Longitudinal stick Sstick, long 
potentiometer 


Left and right 


Throttle position 6+ و‎ бт, spring loaded, modified selsyn 
р potentiometers 
Airspeed Vo W. O. Leonard Model 50-2100 


Altitude h W. O. Leonard Model 50-2000-9 





TABLE 11 
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RECORDED FLIGHT DATA QUANTITIES, 


TELEMETERING CHANNELS UTILIZED 


UANT IT Y 
Roll rate, cg 

Pitch rate, cg 
Yaw rate, cg 
Roll rate, seat. 
Pitch rate, seat. 
Yaw rate, seat. 
Pitch rate, vertical gyro 
Roll rate, vertical gyro 
Lateral acceleration 
Normal acceleration 
Rudder deflection 
Left ai leron deflection 
Left spoiler deflection 
Stabi lator deflection 
Rudder pedal deflection 
Lateral stick deflection 
Longi tudinal stick deflection 
Left throttle 
Right throttle 
Airspeed 


Altitude 


CHANNELS UTILIZED 
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TABLE 111 
WIND CONDITIONS DURING CARRIER TEST APPROACHES 


(Headwind and Crosswind Components) 


RELATIVE RELATIVE 
APPROACH NUMBER HEADWIND (KNOTS) CROSSWIND (KNOTS) 
| ጋዕ 0 
2 26 0 
5 56 0 
4 36 0 
5 36 0 
6 36 0 
7 56 0 
8 19 0 
9 19 ሀ 
10 19 0 
|| 29 ሀ 
[2 26 0 
E 37 0 
14 3| 0 
[5 28 0 
l6 26 0 
|7 25 0 
18 25 0 
19 25 0 
20 25 0 
21 25 0 
22 25 0 
25 25 0 
24 26 0 
25 24 7 1/2 right 
26 28 5 right 
27 24 5 right 
28 27 5 right 
29 55 0 
30 34 6 right 
51 54 5 right 
32 54 9 right 
55 30 5 right 


| 
i 
i 
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TABLE IV 
RMS VALUES OF TURBULENCE QUANTITIES 


ад 9g Bg Qg ў 
(degrees) (deg/ sec) (degrees) (dec/sec) (dedi c) 


APPROACH 
| 2.21 3,06 |,94 2,88 3,39 
2 2.21 2, 78 1.69 - - 
3 2.39 2.89 [207 2,77 4.18 
4 2221 22፡1 - 2.89 3.86 
5 2.16 2072 |, 75 3,19 5.55 
6 2:15 2.77 1.66 5. 27 3,75 
- 2。16 3,17 1.79 5,27 5.73 
9 
1 DATA UNUSABLE 
2 
13 
14 1.64 - 1.86 2.70 5,76 
15 - 3.60 „05 2.71 3.66 
16 1.66 - - 5,54 - 
17 1.63 2.69 ፦ - - 
18 | 5 51. | .97 - в || 
19 |, 42 2,99 |,97 3.03 3.66 
20 - - |, 56 5,54 4.25 
21 - 2.97 ー 2.95 5,54 
22 72152 2.97 1.44 253 4.24 
23 - 3.17 - 3.22 3,66 
24 225] 3.14 1.61 3.51 3.89 
25 2.41 2.89 .57 2.95 5, 48 
26 2,09 3,34 - 3,04 5,78 
27 - 3.16 1.77 - 5,21 
28 = = 2 = = 
29 2.51 = |,67 2,84 - 
30 2,25 Y |, 98 5.14 5,76 
51 2.09 5.27 1.91 2.86 5,89 
52 2,45 ӘЛІ? 2.10 2.99 5.67 
55 2.20 5,27 2.10 5,04 3.83 
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